LPT-ORSAY 02-100 
FAMN SE-02/02 
CPHT RR 089.1202 
UHU-FT/02-24 



(N 
O 
O 
(N 

o 

Q 

^' The strong coupling constant at 

small momentum as an instanton 
(n: detector. 

(N 
^ ■ 

Csl ■ Ph. Boucaud'', F. De Soto% A. Le Yaouanc'', J. P. Leroy'', J. Micheli'', 
^ ! H. Moutarde'^, O. Pene°, J. Rodn'guez-Quintero'^ 

^1 Laboratoire de Physique Theorique Q 

Universite de Paris XI, Batiment 210, 91405 Orsay Cedex, France 
^ . ^Dpto. de Fisica Atomica, Molecular y Nuclear 

Universidad de Sevilla, Apdo. 1065, 41080 Sevilla, Spain 
I Centre de Physique Theorique Ecole Polytechnique, 91128 Palaiseau Cedex, 

I France 

'^Dpto. de Fisica Aplicada e Ingenieria electrica 
E.P.S. La Rabida, Universidad de Huelva, 21819 Palos de la fra., Spain 



Abstract 



We present a study of omomCp) at small p computed from the lattice. 
It shows a dramatic oc law which can be understood within an instanton 
liquid model. In this framework the prefactor gives a direct measure of the 
instanton density in thermalised configurations. A preliminary result for 
this density is 5.27(4) fm"'^. 
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1 Introduction 



In a series of lattice studies the gluon propagator in QCD has been com- 

puted at large momenta, and it was shown that its behavior was compatible with 
the perturbative expectation provided a rather large l/p^ correction was consid- 
ered. In an OPE approach this correction has been shown f^, ^ to stem from an 

gluon condensate which has not to vanish since the calculations are performed 
in the Landau gauge. 

In the deep infrared (IR) region, where the perturbative approach is com- 
pletely meaningless, the present knowledge of the coupling constant is not so 
clear, in spite of all the effort that has been dedicated for years (See and 
references therein). 

Lattice calculations provide a unique laboratory to obtain the behavior of the 
coupling at scales below the hadronization scale, where a complete understand- 
ing of the non-perturbative coupling would be a very important step towards the 
comprehension of hadronization. In this framework, the non-perturbative cou- 
pling has been computed from different vertices, the quark-gluon vertex 0, the 
ghost-gluon vertex p, the three-gluon vertex [y, 0, different propagators 0, 0, 
etc. 

Instantons |n|] have been proposed, via the instanton hquid picture, to de- 
scribe a number of non-perturbative phenomena in QCD, and particularly to 
explain the low part of the Dirac operator spectrum and hence the genesis of the 
chiral Goldstone boson, see |]12|-||15[ to mention only a few papers. 



Instanton studies on the lattice have used essentially the cooled gauge con- 



figurations, [T^-[^, which allow to see instanton-like structures. These results 



have been used to introduce a new definition of the strong coupling constant in 



the IR 121 



This cooling method has been criticized p3[ as creating a distortion on the 
original thermalised gauge configuration, and it was claimed that a direct study of 
the local chirality on a thermalized gauge configuration contradicted the instanton 
liquid picture. However other authors, P^ , , concluded from a similar analysis 
that the dominance of instantons on topological charge fluctuations is not ruled 
out by local chirality measurements. 

In this letter we study another observable, namely the strong coupling constant 
ctMOM the deep IR region on thermalised gauge configurations and we show that 
it strongly supports the instanton liquid picture. 

Simultaneously we think that we provide a very simple and appealing under- 
standing of the IR behavior of omom- This tends to confirm the claim presented 
in 1^ that an instanton liquid might explain the < > condensate observed 
via power corrections to the perturbative behavior of the gluon propagator and 
«MOM in the large momentum regime. 

We will recall the lattice definition of omom, derive the behavior of omom 
in an instanton liquid and compare the latter with numerical results in the low 
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momentum region. We then use cooled gauge configurations to compare, as a 
test, the instanton density derived from aMOM to that which is directly observed 
from shape recognition. We then conclude. 



2 Instanton background effect on aMOM 
2.1 Non-perturbative definition of aMOM 

Let us recall shortly the non-perturbative MOM definition of ^^(p^) |^ in 
Landau gauge. We consider the three-gluon Green function G^^^'^^^"i^^^^{pi, p2, Ps) 
at the symmetric point, Pi = pi = = Z^^- 

The tree-diagram three-gluon vertex is given by Qs T*^^*^ with T*'''^^ defined by 



^m'^Ms = ['^MiM^(Pi -P2);.^ + cycl. perm.] [] ('^m^a 

i=l,3 \ 



Pi fi'-Pi /li 



Pi 
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The three-gluon Green function may be expanded on a basis of tensors. We 
are interested in the scalar function /x^, /i^) which multiplies T**"^*^. It is 

obtained by the following contraction 



,lr,« , (Pl - P-l)n (P2 ' Ps),., (P3 - Pi )m 



rptree , 

I In llnllo 



(2) 



The Euclidean two point Green function in momentum space writes in the 
Landau Gauge: 

G^'Z:. (P^ -P) = G^'' (P'X^^ (^M..2 - ^) (3) 

where oi, 02 are the color indices ranging from 1 to 8. 

Then the renormalised coupling constant is given by [^] 

^"^^^ ' G'(2)(p2)G'(2)(p2)G'(2)(p2) W 

where 

Z3(/i2) = G(2)(/iV (5) 
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2.2 Solution in an instanton liquid 

Let us consider an instanton liquid picture [|12|-||15| i.e. a gauge field given by 

4'^'^(^) = E ^^^(^^ - 4) P7'P i^-^^P^) ' (6) 

i \ Pi / 

where {pi) are the center (radius) of the instantons, r/°^ is known as 't Hooft 
symbol, are color rotations embedding the canonical SU{2) instanton into 
the SU{3) gauge group, a = 1, 3 (a = 1, 8) is an SU{2) {SU{3)) color index, and 
the sum is extended over instantons and anti-instantons. 

In order to take into account instanton deformation resulting from their inter- 
action we do not assume that the radial function P{u) is equal to 2/{v?{v? + 1)) 
as for 't Hooft-Polyakov's instanton. 

The field's Fourier transform is 



^'^'^(P) = ^ E Rl?^%.e^'-'' pf-I{PPi) (7) 

i -P 

for p and ^[/■'"(O) = 0. The function I{ppi) is given by 

4.77-2 „oo 

I{s) = / z^dzJ2{sz)P{z) ; (8) 

s Jo 

J2 being the second order Bessel J function. 
From eqs. (|^), (0) and the relation 

G^'Op) = ^ < 4 (P)4(-P) > (9) 
where the average < • ■ ■ > stands here for the average over instantons, we get 

G^'\p') = l<p'i{ppr> (10) 

n being the instanton density and p being the instanton radius. When computing 
the result in eq. ([T0|) we have used the relations concerning the rj tensors from 
the appendix A in ref. and the orthogonality relation 

Y: Ri^Rt = < E Ri-)R(^)e''^''-'^^ >.^.^ (11) 

a a 

For brevity we skip the derivation of the first of these equations. The second one 
is not exact but expresses the hypothesis ^ that the color rotations are at random 

^We assume a non-negligible instanton density. 

•^In ref. using variational methods it is shown that the color orientation of different 

instantons seems to be weakly correlated. 
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for different instantons as well as the phases e^^'-^*^^^-'. The cross-products are 
thus reasonably assumed to average to zero. 
Similarly, 



G^'\p\p\p') = ^<p'l{ppf> (12) 
48 p 



where we have used again ref. ||T2[ and 

^i?^;^i?gi?g/,,, = e"^^. (13) 



We skip again the derivation of this equation and again we neglect crossed terms 
for the same reason as in (pi]). The strong coupling constant is then given by 



aMOM(p) = -r = T~ 
An An 



Gi^){p\p\p'^ 
{G(^){p^))3 



-(p2G(2)(p2))3/2 



1 



IStt 



(14) 



where we have assumed, for simplicity, that all the instantons have equal radii; 
this will be discussed later. The very remarkable feature of this result is that it 
does not depend on the shape of the instanton-like structures i.e. on the function 
/ (that has all the information about the profile function P{u)), neither on the 
scale p appearing in i.e. on the instanton radius. It only depends on the 
instanton density n. 



3 Lattice three gluon vertex. 

In several lattice works, the scheme outlined in ^TT] has been used to compute 
the running coupling constant in a quite wide range of lattice volumes and spac- 
ings. The region of high momenta has been precisely described according to the 
perturbative running plus non-perturbative corrections in an Operator Product 
Expansion framework (Figure 0(b)). 

The smaller momenta, on the contrary, are not yet theoretically understood. 
The aim of this paper is to try an instanton interpretation at very low momentum, 
as far as possible from the perturbative regime. We will therefore make a fit of 
the deep IR region combining all points from different lattice settings in order 
to achieve enough statistics. The use of different lattice settings, with varying 
statistics, makes it difficult to estimate the errors with the usual jacknife method. 
Therefore, calling x^jj^ the minimum x^, we estimate the errors by assuming that 
one standard deviation is reached when = Xmin + 1 while, of course, the central 
value corresponds to = Xmin- 

A quick inspection of the points in fig. [I|(a) shows that even for very small 
momenta the scaling is very good, i.e. the simulations with different lattice 
spacings agree strikingly. We have used unusually low values of (3 i.e. unusually 
large lattice spacings - and hence large volumes - in order to reach small momenta. 
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The quality of the scahng for these small momenta as well as for larger ones 
makes us confident that we do not introduce a sizable bias with these small /5's. 
We take the ratios between lattice spacings for different /5's from ref. |27| and 
a-\(3 = 6.0) = 1.97 GeV. 

If we fit the low tail of the coupling (up to momenta ~ 0.8 GeV) by a power 
law, a power 3.84(8) is obtained, in good agreement with the expected power-four 
behavior. 

Fixing now the power to be four, the density resulting from formula (0) 
is n = 5.27(4)fm~ , with x^/d.o.f. = 3.8 (Figure ||(a)). This is in the right 



ballpark since several arguments point toward a few fm 's. The rather 



large x^/d.o.f. = 3.8 is due to two points which have unusually small errors, 
maybe because of the small statistics used in this preliminary study. Furthermore 
an exact power 4 should not be taken too seriously since the instanton radii 
distribution can presumably distort this power law. Indeed the dispersion of the 
't Hooft instanton radii leads to an effective power smaller than four at small 
momenta, but the deviation from four is small. 

So within this approach, we are able to compute the instanton density directly 
from the thermalised lattice, and we obtain a result which is not biased by a cooling 
procedure. 





(a) 



(b) 



Figure 1: (b) Symmetric MOM coupling constant for different lattice settings and fits 
to perturbative expression plus power corrections in the high momenta region and to 



expression (14) discussed in the text for small momenta, (a) Region of small momenta 
is zoomed. 



An important question arises here. How can we be sure that what we see 
are really instanton-like objects ? We have already noticed that the result in 
(|I4D is obtained for any radial profile of the semi-classical structures considered, 
provided the tensorial structure is that of eq. (|^). But it is easy to see that any 
semi-classical field configuration different from eq. (P) would also produce a p^ 
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law but with a different prefactor. Therefore we might as well have seen other 
structures than instantons and thus our estimate of the density, which relies on 
the prefactor as compared to the one-instanton prediction, could be wrong. 

In order to check our interpretation we appeal to cooled configurations. We 
insist that we do not use the controversial cooled configurations to infer the 
properties of the thermalised ones, but only to check our new method to estimate 
the instanton density against the instanton shape recognition (ISR) method which 
can only be applied after cooling. We cool the configurations according to the 
method described in p and compute the "strong coupling constant" according 
to the scheme described section The first striking result is that after cooling 
the coupling constant grows in the whole energy range (Figure It reaches 
exceedingly large values at large momenta, the reason of which being the oc 
law and the large prefactor due to a small instanton density, see eq. (|T^) . 

We now make a fit to a p'^ law for auouip) on lattice configurations which have 
been submitted to a large number of cooling sweeps (200) so that UV fluctuations 
should have disappeared leaving only semiclassical structures. We do not take in 
the fit the few smallest momenta for reasons which will be discussed soon. We 
compare the instanton density extracted via formula (0) applied to the oc p^ 
fit of auouip) to the instanton density coming from the geometrical "instanton 
shape recognition" (ISR) method described in p. We obtain the qualitative 
agreement shown in table |1|. Notice that these values for the instanton densities 
have nothing to do with the latter density in the thermalised configuration Q 
What is relevant is the fair agreement between both methods to estimate the 
instanton density. 



L 


P 


n(ISR) 


n(a) 


24 


5.6 


0.009(6) 


0.019(1) 


24 


5.8 


0.048(23) 


0.090(1) 


24 


6.0 


0.115(6) 


0.133(2) 


32 


6.0 


0.145(16) 


0.197(7) 



Table 1: Comparison between the instanton density, n (in fm~^), obtained through 
the Instanton Shape Recognition (ISR) method and the density deduced from the 
fit of the coupUng constant after 200 cooUng sweeps. In both cases errors are only 
statistical. 



We expect differences in table |1| between both estimates, mainly at small 
(5, because the ISR method does not recognize small instantons in lattice units, 
which induces for ISR a systematic underestimate of the density p8[ . The general 



^For simplicity we call "coupling constant" this quantity computed according to the defini- 
tions in section 2.1 although being aware that in cooled configurations this denomination is not 



really appropriate. 

^For simplicity we have used the lattice spacing of the thermalised configurations. 
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tendency in table |1] supports this argument, as both estimates of the density are 
closer for large values of beta f\. 

In figure ^(a) it can be seen that at low momentum the points for auoM 
are below the fit. The log- log plot of the same quantity for different cooling 
sweeps in figure ^(b) confirms that the behavior at small momenta differs slightly 
from the large momenta one. In fact the detailed dependence of the slopes in 
the momentum and the number of cooling sweeps seems to be involved but it is 
striking that the slope stays always in the range three to five, i.e. close to the 
expected slope four. 

Anyhow, why do the thermal configuration seem to agree at small momentum 
with the oc p"^ law, fig. |l|(a), while the cooled configurations seem unexpectedly to 
agree with the same law only at larger momentum ? One might argue that at large 



distance instanton deformation |28| as well as instanton (anti-)instanton repulsion 
(attraction) might have been generated by the cooling itself []. This would explain 
why the law observed at small momentum in the thermalised configurations 
tends to be distorted with cooling. Honestly this argument has to be submitted 
to closer scrutiny and for now we consider the detailed understanding of these 
slopes as an open question, but we would like to stress that anyhow the observed 
behaviour is never far from oc p^. 

4 Conclusion and discussion 



We have found that the lattice simulations indicate good scaling of aMomip) 
when the lattice spacing is varied, even at rather small (3 when p is small. 



• We find ^ 

auouip) ^ -^ri-y, for p < O.SGeV (15) 

as expected from an instanton liquid picture and eq. ([l^ . 

• The fitted density 

n = 5.27(4)fm-^ (16) 
is in fair agreement with expectations. 

• We have checked that the calculation of the instanton density from eq. (0) 
via a oc fit of q;mom(p) is comparable to a direct counting of instantons 
recognised from their shape. We take this as a convincing evidence that 
this new aMOM-method to measure instanton density is reliable. 



^As the lattice spacing is increased an increasingly large number of instantons are missed 
by the ISR method because they become too small in lattice units. 

^During the cooling instanton anti-instanton pair annihilation occurs which confirms that 
correlations are generated by the cooling. 
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a b 

Figure 2: (a) Coupling constant in a cooled lattice {L = 24, /3 = 6.0) after 200 cooling 
sweeps. The solid line corresponds to the fit discussed in the text. The horizontal 
axis is given in GeV, assuming for simplicity the lattice spacing of the thermalised 
configurations, = 1.97 GeV. (b) Coupling constant in a cooled lattice (L = 24, /3 = 
5.6) after several cooling sweeps. The horizontal axis is the momentum assuming 
the same lattice spacing, = 0.83 GeV, which is the value for the thermalised 
configurations. This is a log-log plot which exhibits better the power law. There seems 
to be three regimes. 



This makes it highly plausible that this aMomip) "i^ C''^ ^ff^ct of a liquid 
of instantons, i.e. that such a liquid of instantons indeed exists in the ther- 
malised configurations and that the quantum fluctuations do not affect signifi- 
cantly auiouip) for p < 0.8 GeV. 

These results open a Pandora box of new questions: can this simple expla- 
nation also apply to other definitions of a.^, to other Green functions (the gluon 
propagator in particular) ? How is this interpretation related to Schwinger-Dyson 
or renormalisation group deduced small momentum behavior of these quantities ? 

A look at figure ^ (b) shows a nice theoretical understanding (solid lines) 
of aMouip) both in the large (perturbative QCD + OPE) and small (instanton 
liquid picture) momentum regimes. How to understand better the transition 
between these two regimes ? 
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